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Graphene’s linear dispersion relation and the attendant implications for bipolar electronics applications have 
motivated a range of experimental efforts aimed at producing p-n junctions in graphene. Here we report elec¬ 
trical transport measurements of graphene p-n junctions formed via simple modifications to a PbZro.2Tio.8O3 
substrate, combined with a self-assembled layer of ambient environmental dopants. We show that the sub¬ 
strate configuration controls the local doping region, and that the p-n junction behavior can be controlled 
with a single gate. Finally, we show that the ferroelectric substrate induces a hysteresis in the environmental 
doping which can be utilized to activate and deactivate the doping, yielding an ‘on-demand’ p-n junction in 
graphene controlled by a single, universal backgate. 
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Graphene is a subject of intense research interest due 
to the enormous potential of its electronic and mechani¬ 
cal propertied. In particular, p-n junctions in graphene 
have great potential for both fundamental research and 
commercial application^ and have been utilized to study 
the quantum Hall effecfP®! and Klein tunnelin^^ as well 
as to fabricate flexible transistord. Previous work on p- 
n junctions in graphene employed multiple electrostatic 
gated^^HI^, charge transfer from the controlled depo¬ 
sition of chemical adsorbateJ^^®^, high current-induced 
charging of trap states in the substratd^, or periodically 
poled ferroelectric substrate^. 

In this Letter we report the fabrication of p-n junctions 
in graphene deposited on a uniformly poled ferroelectric 
substrate. The local doping in our devices is accom¬ 
plished by combining simple modifications to a lead zir¬ 
conium titanate (PbZro.2Tio.8O3) substrate - the evap¬ 
oration of thin Si02 films in some regions - with a self- 
assembled layer of environmental dopants. We find that 
the PbZro.2Tio.8O3 substrate modulates the doping effect 
of the adsorbed dopants: devices are exposed to ambient 
conditions after fabrication whereupon experimental ob¬ 
servations confirm both the presence of adsorbed dopants 
(likely primarily H2O) and their enhanced doping effect 
on the PbZro.2Tio.8O3 relative to the Si02. Further¬ 
more, we demonstrate that the PbZro.2Tio.8O3 substrate 
induces a hysteresis in the environmental doping which 
can be used to activate and deactivate the doping via the 
application of large gate voltages. We employ this effect 
to create p-n junctions which can be reversibly transi¬ 
tioned between p-n junction and uniformly conducting 
configurations. 

Devices consist of graphene micro-ribbons deposited 
on substrates which are partially covered by a thin layer 
of evaporated Si02, and contacted in a four-point geom¬ 


etry, as illustrated in Figures [^-c. An SEM micrograph 
of a typical device is shown in the inset of Figure p. 
The devices are fabricated using standard lithography 
and deposition techniques on thin-fllm ferroelectric sub¬ 
strates. For the ferroelectric substrates, 120 nm thick 
(OOl)-oriented lead zirconium titanate (PbZro.2Tio.8O3) 
films are prepared by pulsed-laser deposition (PLD) on a 
strontium titanate (SrTi03) substrate coated with 20 nm 
of stronti um rn thenate (SrRu03), following established 
procednre j^^ l ^^ i For each substrate, an 80 nm-thick layer 
of Si02 is evaporated in small rectangular regions, with 
region widths ranging from 0.5 /im to 3 /im, as illustrated 
in Figs. P-c. GVD graphene is then transferred using 
standard wet transfer techniqueP^, and patterned into 
ribbons spanning the deposited Si02 using photolithog¬ 
raphy and reactive ion etching. The graphene channel is 
6 /im by 4 /im measured from the inner contacts. Gon- 
trol devices which span regions with no evaporated Si02 
are also fabricated. Finally, Gr/Au (3 nm /20 nm) leads 
are deposited in a four-point measurement configuration. 
Raman spectroscopy is used to confirm the monolayer 
character and high quality of the graphene after fabrica¬ 
tion; a representative spectra is shown in Figure [^. 

Transport measurements in air are performed using 
two Keithley 2400 SourceMeters. Measurements in vac¬ 
uum are performed using an Agilent 4156 G Semiconduc¬ 
tor Parameter Analyzer. In both cases source-drain cur¬ 
rent is measured with a constant source-drain bias of 5 
mV while the voltage applied to the backgate is swept. 
Gate leakage distorts transport results for gate voltages 
more positive than 1.5 V or more negative than -IV, 
so gate voltages are limited to this range. Gate voltage 
sweep rates range from 10 mV/s to 100 mV/s; the data 
presented here are from sweeps at 100 mV/s. Slower 
sweep rates yield qualitatively similar results. 
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FIG. 1. (a-c) Schematic illustration of the device geom¬ 

etry. The graphene (grey) spans the regions of evaporated 
Si 02 (purple) on the bare substrate (teal) and is contacted by 
Cr/Au electrodes (yellow). The universal backgate is shown 
in red. The evaporated Si02 region width ranges from 0.5 /xm 
to 3 /xm in the measured devices. The graphene channel is 6 
/xm by 4 /xm measured from the inner contacts, (d) Raman 
spectra of graphene from a device fabricated on a thermal 
Si02 substrate. Inset: a false-color SEM micrograph of a 
representative device. The scale bar is 2 /xm. 


Figure]^ shows room temperature Isd vs Vgate curves 
for devices having different widths of evaporated Si 02 
on a PbZro. 2 Tio. 8 O 3 substrate. For the data shown here, 
the fraction of the graphene channel screened by evap¬ 
orated Si02 ranges from 0% to 50% (corresponding to 
Si02 widths of 0 to 3 microns). Two features of the 
data are immediately apparent: first, for devices which 
span an evaporated Si 02 region, the characteristic con¬ 
ductance minimum typically observed in graphene at the 
Dirac point is split into two distinct minima, one at the 
original Dirac point location and a second shifted to the 
right. This is apparent in the top and bottom curves of 
the Figure: the bottom curve, corresponding to a device 
with 0% screening, displays a single minimum, while the 
top curve, corresponding to a device with 50% screen¬ 
ing, displays two pronounced minima. Second, the width 
of the evaporated Si 02 region determines which of the 
two minima has a smaller absolute value. As the screen¬ 
ing fraction is increased, ‘weight’ is transferred from the 
minimum at the original Dirac point location to the sec¬ 
ondary minimum, and the depths of the two minima vary 
accordingly. We attribute both of these effects to the 
presence of two different doping regions in the graphene, 
defined by the PbZro. 2 Tio. 8 O 3 substrate and the evapo¬ 
rated Si 02 . 

The data can be understood by considering that as the 
gate voltage is swept from negative to positive, the Fermi 
level passes through the charge-neutrality point (GNP) of 
each graphene region separately. Taking the conductance 
to be linear with carrier densit}E31 a oc kpir) oc n and the 
carrier density to be linear with the thermally smeared 


energy difference between the Fermi level and the GNpl^, 
we model the conductance in the vicinity of the GNP as: 

oc n oc 1 — 1 2 c ^ where the constant 

fi accounts for the extrinsic doping introduced by the 
fabrication process, J G {0,1} describes the substrate- 
dependent doping, and e accounts for the non-vanishing 
carrier density at the GNP. Assuming diffusive transport 
in the graphene, the relative weight of each separately 
doped region, and therefore the relative magnitude of 
the measured conductance minima, is determined by the 
fraction of the graphene channel which is screened: 

Gd OC [pscr. X (Pcf-scr.) T Pnon-scr. X -non-scr.)] 

This is simulated in Figure]^ which shows Isd vs Vgate 
curves for screening fractions ranging from 0% to 75% 
and is in excellent agreement with our experimental data. 
We note that the simulations agree with our data for 
/X > 0 and J > 0, which is consistent with the extrinsic p- 
doping typically observed in graphene devices fabricated 
on PbZro.sTio.gO#^. 

The substrate-selectivity of the doping in our de¬ 
vices suggests that the ambient dopants are polar H 2 O 
molecules. Polar surface adsorbates have been shown to 



FIG. 2. (a) Offset Isd vs Vgate curves for devices hav¬ 

ing different widths of evaporated Si02 on a PbZro. 2 Tio. 8 O 3 
substrate. The evaporated Si02 serves to locally screen the 
graphene from any substrate-dependent effects; from bottom 
to top the curves correspond to graphene channels which are 
0, 17, 25, 42, and 50 percent screened. The devices hav¬ 
ing evaporated Si02 show split Dirac points, and the relative 
dominance of the left and right Dirac points can be tuned 
by varying the evaporated Si02 width, (b) Simulations of 
transport across a graphene device having two different lo¬ 
cally doped regions, as a function of applied gate voltage. 
The simulations show both Dirac point splitting and varia¬ 
tions in left vs. right Dirac point dominance as a function 
of screening region width, consistent with the experimental 
data. From bottom to top the curves correspond to graphene 
channels which are 0, 10, 25, 50, and 75 percent screened. 
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dramatically affect the electronic pr operties of complex 
oxide systemas well as graphen ^^ l ^^ l ^^ i The H2O 
doping is substrate-selective because of the unique ferro- 
elect ric na ture of the PbZro.2Tio.8O3 substrate. Previous 
worlPSEB Pas established the importance of the orbital 
structure of the adsorbate in determining the most en¬ 
ergetically favorable orientation. For standard graphene 
devices on Si02 substrates, these energetics favor a uni¬ 
form polarization throughout the range of applicable gate 
voltages. In the devices described here, it is likely that 
the remnant polarization of the ferroelectric substrate 
sufficiently alters the orbital structure of the adsorbed 
H2O molecules to destroy this stability, and thus create 
less-polarized regions. This hypothesis is in agreement 
with data collected from similar devices fabricated on 
non-ferroelectric substrates; such devices show none of 
the characteristic Dirac point splitting associated with 
p-n junctions. 

The gate-voltage dependence of the H2O polarization 
configurations on PbZro.2Tio.8O3 vs Si02 leads to hys¬ 
teresis in the devices. This can be seen in figure [^, 
which shows Isd vs Vgate curves for forward and reverse 
gate sweeps performed on the same devices as measured 
in Fig. A pronounced hysteresis between forward 

and reverse gate sweeps is apparent. As in Fig. 1 ^, de¬ 
vices spanning a region of evaporated Si02 display two 
distinct minima during forward sweeps, while a control 
device having no Si 02 displays a single minimum. How¬ 
ever, all devices display a single minimum during reverse 
gate sweeps. We note that the gate voltages applied here 
remain below the coercive voltage of the PbZro.2Tio.8O3 
film, therefore ferroelectric switching is not a possible 
cause of the observed hysteresis. 

In order to understand how the hysteresis is related 
to the different H2O polarization configurations on the 
PbZro.2Tio.8O3 as compared to the evaporated Si02, it is 
instructive to consider the gate voltages at which the var¬ 
ious minima appear. For example, for the 17 % screened 
curve in Fig. the red arrows point out two min¬ 
ima on the forward sweep (at 0.6 V and I.IV) and one 
minimum on the reverse sweep (at 0.9 V). These can 
be compared to the position of the Dirac point in vac¬ 
uum at 0.6 V (see Fig. [^. The minimum on the for¬ 
ward sweep at 0.6 V occurs at the same gate voltage 
as the vacuum Dirac point, implying that it corresponds 
to a region of the graphene without a net polarization 
in the adsorbed H2O. The remaining minima occur at 
voltages larger than the Dirac point ( 0.9 V and 1.1 V) 
and thus correspond to regions of the graphene on which 
the adsorbed H2O is polarized and produces p-doping. 
Polarized H2O typically produces p-doping in graphene, 
though the preci se mechanism is the subject of continu¬ 
ing researWe identify the forward-sweep min¬ 
imum at 0.6 V as corresponding to graphene on the non- 
screened PbZro.2Tio.8O3 region. This is supported by 
the fact that all devices demonstrate a minimum at 0.6 
V, independent of different Si02 screening fractions. In 
contrast, the minimum indicated by the rightmost arrow 



FIG. 3. (a) Offset Isd vs Vgate curves for forward (light) 

and reverse (dark) gate voltage sweeps. In all non-control 
devices, forward sweeps show split Dirac points, but reverse 
sweeps show only a single Dirac point. Representative min¬ 
ima locations are indicated by the vertical arrows, (b) H 2 O 
polarization by device region, as prepared by positive and neg¬ 
ative gate voltages. For negative gate voltages, H 2 O on the 
PbZro. 2 Tio. 8 O 3 region is unpolarized, and H 2 O on the Si 02 
region is polarized. For positive gate voltages both regions 
have a net polarization, (c) Simulated carrier density vs. ap¬ 
plied gate voltage. The onset of H 2 O polarization (indicated 
by grey shading) decreases the slope of the carrier density 
vs. gate voltage curve for graphene on PbZro. 2 Tio.8 O3. Inset: 
simulated conductance vs. applied gate voltage; the slope of 
carrier density vs. gate voltage at the CNP determines the 
width of the conductance minimum, (d) A schematic illustra¬ 
tion of the H 2 O doping hysteresis with applied gate voltage; 
arrows indicate the direction of the gate voltage sweep. The 
application of a large negative gate voltage destroys the H 2 O 
polarization while a large positive gate voltage establishes a 
net polarization. 

(e.g. at 1.1 V for the 17 % screened device) corresponds 
to graphene on the evaporated Si02 region, as evidenced 
by its evolution with increasing Si02 screening fraction. 

We conclude that the application of a negative gate 
voltage destroys the net polarization of adsorbed H2O on 
PbZro.2Tio.8O3, but preserves a net polarization on the 
Si02-screened regions. This creates different local doping 
levels and thus a p-n junction. A positive gate voltage 
establishes a net polarization in both regions, and thus a 
uniform channel with no p-n junction, as depicted in Fig¬ 
ure Dd. This interpretation is further corroborated by the 
different widths of the forward and reverse minima for the 
control ( 0 % screened) device, as evident in Fig. [^. This 
difference can be explained by considering that conduc¬ 
tance is linear with carrier densit}^^, so the width of the 
conductance minimum at the Dirac point is determined 
by the slope of the carrier density vs. gate voltage curve. 
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Typically the slope is constant, determined by the gate 
capacitance. However, for our devices the onset of H 2 O 
dipole doping introduces a nonlinearity in the regime 
where the adsorbed H 2 O transitions from unpolarized to 
polarized; this is shown schematically in Fig. li- The po¬ 
larized H 2 O in our devices p-dopes the graphene, so the 
onset of its doping contribution temporarily reduces the 
slope of the carrier density vs. gate voltage curve, broad¬ 
ening the conductance minimum. During reverse sweeps, 
the transition from polarized to unpolarized H 2 O occurs 
far from the GNP, so the width of the conductance min¬ 
imum is unaffected. The hysteresis in H 2 O polarization 
is illustrated schematically in Figure [^. Experimen¬ 
tally, the polarization hysteresis displays a dependence 
on both the magnitude of the applied gate voltage and 
the duration of its application, which prevents an exact 
determination of the gate voltages required to establish 
or destroy the H 2 O polarization. 

The hysteresis of the H 2 O polarization on 
PbZro. 2 Tio. 8 O 3 substrates adds an ‘on/off’ switch¬ 
ing element to the p-n behavior. Specifically, we can 
selectively transition the device into and out of the 
p-n junction configuration through the application of 
large positive and negative gate voltages. The initial 
application of a large positive gate voltages establishes 
a uniform polarization across the device, yielding a 
unipolar conducting channel, while a large negative 
gate voltage destabilizes the polarization on regions 
supported by PbZro. 2 Tio. 8 O 3 . In the latter case the 
different H 2 O polarizations create separate locally doped 
regions, and thus a p-n junction. This ‘on/off’ switching 
is different than the standard gate induced switching 
observed in p-n junctions, for example from p-n to p^-p. 
By comparison, in our devices the same applied gate 
voltage can generate either a p-n junction or a uniformly 
doped channel, depending on the H 2 O polarization 
condition. 

The ferroelectric nature of the PbZro. 2 Tio. 8 O 3 sub¬ 
strate might suggest that the residual electric field from 
the substrate polarization dopes th e regi ons of graphene 
in direct contact with the substrate^^^, but has less ef¬ 
fect in the graphene regions screened by evaporated Si 02 . 
However, this explanation is precluded by several further 
experimental observations. First, the Dirac point split¬ 
ting effect disappears when the devices are measured in 
vacuum. Figured shows Isd vs Vgate curves for the same 
devices measured at 5 x 10“^ Torr but otherwise under 
conditions identical to those of Figure 1^. All devices 
show a single minimum, independent of evaporated Si 02 
width or gate sweep rate. Leaving the devices in ambi¬ 
ent conditions overnight recovers the splitting effect. The 
observed behavior is consistent with an ambient dopant 
mechanism, i.e., the substrate-selective formation of a 
self-assembled layer of dopant molecules. In vacuum, 
dopant molecules desorb from the surface leaving all re¬ 
gions of the graphene identically doped. Leaving the de¬ 
vice in ambient conditions allows the dopant layer to re¬ 
assemble, thereby re-establishing the separately doped 



FIG. 4. Offset Isd vs Vgate curves for devices measured in 
vacuum. The Dirac point splitting effect disappears in vac¬ 
uum, but can be recovered by leaving the device in ambient 
conditions overnight. From bottom to top the curves corre¬ 
spond to graphene channels which are 0, 25, and 50 percent 
screened. 


regions. 

The absence of Dirac point splitting in vacuum mea¬ 
surements also eliminates differences in gate capacitance 
as a dominant source of the splitting effect. The screened 
and non-screened regions of the device have different gate 
thicknesses and dielectric constants, which might suggest 
that the application of the same gate voltage would gen¬ 
erate different doping levels in each region, and hence the 
transport behavior we observe. However, any capacita- 
tive differences between the regions are static, depend¬ 
ing only on the geometry of the device, while the Dirac 
point splitting effect is dynamic, disappearing in vacuum. 
Gapacitance-based explanations are further ruled out by 
the nearly identical Dirac point locations observed in all 
devices under vacuum. In particular, the similarity of 
data under vacuum from the control device (having no 
evaporated Si 02 ), and from the devices which do span an 
evaporated Si 02 region confirms that gate capacitance 
differences between the two regions are not the primary 
cause of the Dirac point splitting effect. 

In summary, we have fabricated a controllable p-n 
junction in graphene on a ferroelectric substrate. We em¬ 
ploy simple substrate modifications to define local doping 
regions, where the doping is accomplished through the 
substrate-selective formation of a self-assembled layer of 
ambient doping molecules. Alternative explanations for 
the local doping effect are ruled out, and the dynamics of 
the ambient doping suggest that it is due to polar H 2 O 
molecules. Finally, the PbZro. 2 Tio. 8 O 3 substrate creates 
a hysteresis in the ambient doping effect which can be 
used to controllably bias the device into and out of a p-n 
junction configuration, using a single, universal backgate. 
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